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Abstract 

We propose a systematic method to extract conformal loop models for rational con- 
formal field theories (CFT). Method is based on defining an ADE model for boundary 
primary operators by using the fusion matrices of these operators as adjacency matrices. 
These loop models respect the conformal boundary conditions. We discuss the loop mod- 
els that can be extracted by this method for minimal CFTs and then we will give dilute 
0(n) loop models on the square lattice as examples for these loop models. We give also 
some proposals for WZW SU(2) models. 
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1 Introduction 

The study of statistical models related to loop models is interesting both from the physical 
and the mathematical point of views. Most of the statistical models studied in physics such 
as the Ising, the q-state Potts model and also complicated vertex models can be represented 
in terms of loops pQ. The loop representation of the spin system is very easy to understand: 
loops correspond to domain walls separating regions of different magnetization. The study of 
critical loop models can be interesting from many point of views: they are good candidates 
for the ground state of topological quantum systems [2], they are also good candidates for 
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the Schramm Loewner evolution (SLE), a method discovered by Schramm [3] to classify 
conformally invariant curves connecting two distinct boundary points in a simply connected 
domain. 

Different applications of conformal loop models are stimulating to do a systematic study 
of these models by CFT. Recently we proposed in [1] a method to extract loop models 
corresponding to a conformal field theory (CFT), the method was based on defining a RSOS 
model for every primary operator by using fusion matrix of the primary operator as an 
adjacency matrix and then extracting the loop model corresponding to domain walls of the 
RSOS model. The weight of the loop model is equal to the quantum dimension of the 
corresponding operator. In this paper we want to follow the same method consistent with 
the conformal boundary operators, since the SLE is a boundary CFT we think that using 
the fusion matrix of boundary operators as an adjacency matrix is more consistent with the 
nature of SLE. Recently a very nice and strong project was initiated by Jacobsen and Saleur 
[5] followed by Dubail, Jacobsen, Saleur [6] to classify all the possible conformal boundary 
loop models. It is based on classifying the possible boundary loop models compatible with 
the boundary conformal field theories. This classification is in close relation with the earlier 
work by Cardy on formulating the modular invariant partition function of 0(n) model on the 
annulus [7j. The results that we get by our method apart from simplicity are all compatible 
with the results in El [7] . 

The paper is organized as follows: In the next section we will introduce the necessary 
ingredients to find the boundary operators and also the fusion matrices corresponding to 
them. In the third section we briefly review the method proposed in [3] and we will also 
generalize it to the graphs with largest eigenvalue bigger than two. The central claim of 
this section is as follows: the loop model extracted with this method is connected with the 
properties of the statistical loop model in the same universality class as the corresponding 
CFT. In the third section we follow explicitly some examples in particular; Ising model, tri- 
critical Ising model, three states Potts model and tri-critical three states Potts model. Then 
we will give the possible loop models, extractable with this method, of minimal CFTs and also 
the lattice models corresponding to these loop models. We will close this section by giving 
some proposals for possible loop models for WZW SU(2) models. Last section contains our 
conclusions with a brief description of the work in progress motivated by these results. 
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2 Boundary conformal field Theory 



To define loop model for a generic minimal CFT consistent with the conformal bound- 
ary we need to first summarize the main important facts about boundary CFT. The most 
important ingredient to classify the boundary conformal operators is the modular invariant 
partition function of the CFT. The classification of modular invariant partition functions of 
SU (2) minimal models are well known and can be related to a pair of simply laced Dynkin 
diagrams (A,G) [8]. The complete classification based on ADE diagrams is 

(A h -i,D(g +2 )/2), 9 even 

(A,G) = l (Ah^Ek), <? = 12 (2.1) 

(A h ^,E 7 ), 5 = 18 

^A^Es), 5 = 30, 

where g and h are the Coxeter numbers of A and G with h,g > 2. The above pair of Dynkin 
diagrams describes bulk modular invariant partition function with some primary operators 
and with the following central charge 

c = 1-6 V , '- . (2.2) 
hg 

Each of the unitary minimal models M(Ah-i,G) with g — h = ±1 can be realized as the 
continuum scaling limit of an integrable two-dimensional lattice model at criticality, with 
heights living on the nodes of the graph G. In particular, the critical series with g — h = 1 
is associated with the A-D-E lattice models [9] and the tri-critical series with g — h = — 1 
is associated with the dilute lattice models |10t lllj . For theories with a diagonal torus 
partition function it is known that there is a conformal boundary condition associated to 
each operator in the theory [12]. The fusion rules of these boundary operators are just given 
by the bulk fusion algebra. It was shown in a series of papers that for SU (2) minimal models 
one can propose a complete set of conformal boundary operators i = (r,a) £ (A,G), where 
r and a are nodes on the Dynkin diagram of A and G respectively with the identification 
(r, a) = (h — r, 7(a)), where 7 is an automorphism acting on the nodes of the graph G. 
This automorphism is identity except for the A, Eq and D dd which is Z2 symmetry of 
Dynkin diagram, symmetries of Dynkin diagrams play an important rule in the forthcoming 
discussion. Following [13] we show the corresponding operators by fa and the independent 
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boundary states by |(r, a)) which is called Cardy states. Cardy states can be written in 
terms of Ishibashi states, i.e. \j)), as follows |(r, a)) = X^cfJi)), where sum is over all 
Ishibashi states. We are interested to the fusion rules of these boundary operators. To give 
a formula for the fusion rules of these operators we need to define some quantities. Let be 
the eigenvectors of the adjacency matrix corresponding to the group G then the graph fusion 
matrices N a with a G G can be defined as follows 



(Na) 



E 

meExp(G) 



a, b, c G G, 



I'm 



(2.3) 



where Exp(G) denotes the set of exponents of G, see table 1. Let's show also the graph 
fusion matrix for Ah-i by N r then following [13] the fusion rules for boundary operators are 

(2.4) 



i 3 e(A,G) 



where (A/^)-^ has the following relation with the graph fusion matrices of A and G 



(2.5) 



For more details about the connection of the boundary operators to bulk counterparts see 

[HE]. 

To calculate the fusion matrices of boundary operators we need also to define a conjugation 



Dynkin Diagram 


Coexter Number(/i) 


Coexter Exponent (m) 


A n 


71 + 1 


l,2,...,n 


D n 


2(n + l) 


1, 3, 2n — 1, n — 1 


E 6 


12 


1,4,5,7,8,11 


E 7 


18 


1,5,7,9,11,13,17 


Es 


30 


1,7,11,13,17, 19,23,29 



Table 1: The Coexter number and the Coexter exponents of Dynkin diagrams. 



operator C(a) = a*, it is the identity except for D± n graphs where the eigenvectors \I/ am 
are complex and conjugation corresponds to the Zi Dynkin diagram automorphism. It then 
follows that N£» b = N^ a . The operator C{a) acts on the right to raise and lower indices in the 
fusion matrices N a = N a C so it is the important ingredient to get the right fusion matrices 
for the boundary operators, in particular for the D^ n graphs, we will give some examples in 
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section 4, in particular we use the above method to get the fusion matrices of the boundary 
operators of Ising model, tri-critical Ising model, 3 state Potts model and tri-critical 3 state 
Potts model. 

3 Loop Models for Boundary operators 



In this section we propose a method to extract some possible loop models for CFTs, the 
method is the same as the method introduced recently in In that reference we showed 
that using the fusion matrix as an adjacency matrix it is possible to associate a 0(n) loop 
model to every primary operator. The method is briefly as follows: The graph of a primary 
operator fa has g vertices where g is the number of primary operators in the theory and 
edges connecting pairs of vertices (j, k) when Mfj = 1. Following [24] one can define a height 
model on the triangular lattice by imposing that the height hj at the site j can take values 
0,1, ... ,5 — 1. Then constraint the heights at neighboring sites according to the incidence 
matrix associated to a given primary field fa: only neighbor heights hj and with (Mi)j = 1 
are admissible. For a consistent definition of loop models on a triangular lattice at least two 
of the heights at the corners of an elementary triangular plaquette should be equal then the 
weights for the elementary plaquette are defined as follows. If the heights of plaquette are 
(c,b,b) with c / b then the weight is x(-£) ' , where S satisfies J2b(N a )bTo = tt^o- It 

S a 

means that the b th element of the eigenvector of M a with eigenvalue -4- is given by If 
the heights are all equal then the weight is 1 except for those with Mi\ 7^ that have weights 
1 or x depending on the particular model considered Q. The next step is to mark triangles 
with unequal heights (c, b, b) drawing a curved segment on the dual honeycomb lattice [23] 
and linking to the center the midpoints of the two edges with different heights (b and c) at 
the extremes (See fig 1). Summing over the admissible values of heights consistent with a 
given loop configuration we find 

qb da 
b a i °0 

where the sum is just over b. We take most of the times I = to get the largest eigenvalue 
of M a to guaranty positive real weights in our height models, however, we will also point to 
other cases. The weight of the loops is given by the largest eigenvalue of the fusion matrix 
1 For more details specially about identical neighbor heights see [3]. 
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Figure 1: A triangular plaquette with c ^ b and the corresponding curve segment on the dual 
honeycomb lattice. 

n a and the partition function of the model is as follows 

Z = 5>^, (3.2) 

where I is the number of bonds in the loop configuration and N is the number of loops. 
Using this method we can correspond to every boundary conformal operator a O(n) loop 
model, since the 0(n) model posses a dilute critical point for n < 2 with x c = , 1 = : 

V 2+V2— n 

see [15] . correspondingly our loop models will have a critical point just for the fields with n a 
smaller than 2. The 0(n) model has another critical regime, the so-called dense phase, for 
x = (x c , oo) corresponds to a different universality class. Mapping to the 0{n) model helps 
us to find the connection with SLE: from coulomb gas arguments we know that, in the dilute 
regime, the loop weight has the following relation with the drift in the SLE equation 

n a = — 2cos( — ). (3-3) 

K 

For the dense phase the above equation is still true if we work in the region 4 < k < 8. Using 
the above equation we can find the properties of the loop model corresponding to a boundary 
conformal operator. The achievement of this method is respecting the Cardy's equation |12j : 
fields in the same sector have the same loop representation. 

Before generalizing the definition to more general graphs we should stress that although 
we started with well defined minimal CFT but the loop model that we extracted is not 
necessarily minimal. The point is that the extracted loop model respects some aspects of the 
corresponding conformal field theory. This is like to say that although the domain walls in 
Ising model at the critical point is the same as the critical 0{n = 1) but the Ising conformal 
field theory does not explain all the aspects of the critical curves. From this point the loop 
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model that one can get by this method from the rational CFT is not perfectly equal to the 
corresponding CFT. 

One can generalize the above idea to the decomposable fusion graphs by the method 
that was explained in [16]. Since the fusion graphs of some operators in minimal models 
are equivalent to the tensor product of two adjacency diagrams one can use this method to 
extract new loop models that can also have configurations with crossing loop segments. The 
general strategy is based on extracting critical loop models with n < 2 for the graphs with 
largest eigenvalue bigger than 2. Some graphs obey simple decomposition, can be written 
as tensor product, but others need to be mapped to simple decomposable graphs by going 
to the ground state adjacency graph [16]. Here we just comment on decomposable graphs 
M = M\ <8>A/i, where Ai and N% are simple ADE diagrams. In these cases we can define two- 
flavor loop model living on the honeycomb lattice independently, one is related to the loop 
model of Mi with weight ni and the other comes from the graph A2 with weight n,2- Fendley 
showed [16] that in this case it is also possible to define consistently interacting loop model 
on the square lattice with partition function Z = Yl, n^n^ 2 ^ , where N\2 are the numbers 
of each kind of loop and C is the number of plaquettes with a resolved potential crossing at 
their center. The critical values of b were calculated in |17] but the critical properties of the 
loops are still unsolved. This is obviously is not the only method to define loop model for 
non-simple graphs, the other method is based on the multi-flavor loop model of [11]. In this 
loop model a curve of flavor i separating two neighboring sites does not necessarily separate 
two sites with different heights, for the definition of the RSOS model in this case and its 
relation to the loop model see |llj . 

In the next section we summarize some simple examples including the most familiar 
minimal conformal models such as Ising, tri-critical Ising, 3-state Potts model and tri-critical 
3-state Potts model. The main point is to take the fusion graphs as adjacency graphs in the 
consistent way and to extract some loop models. These loop models are not equivalent to 
the corresponding conformal field theory but still carry some aspects of the underlying field 
theory in the consistent way, in particular the critical properties of these loop models are in 
close connection with the corresponding conformal field theory. 

In this paper some distinctions are crucial. We have some minimal conformal field theories 
with well defined fusion matrices and modular invariant partition functions, one example is 
Ising conformal field theory. There are some statistical models such as spin models, RSOS 
models which at the critical point can be describe partially by the minimal CFT, so the 
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Figure 2: Graphs of fusion matrices of boundary primary operators in Ising model, from the 
left to the right the fusion graphs of 1, e and a. The graph of the operator e is A2 and the 
graph of a is ^3 

Ising CFT is different from the statistical Ising model. We prefer also to distinguish between 
for example dilute ADE models and dilute 0(n) loop model. They can be mapped to each 
other and have the same phase transitions but since the fundamental objects in one side 
is local and in the other one is non-local this distinction is useful. There are lots of work 
done on connecting these two models, minimal conformal field theories and statistical models 
counterparts, using integrability methods and our argument hardly has something new to 
say from this point of view. Finally we are defining another statistical model by using the 
fusion matrices of primary operators of conformal field theory which most of the times is in 
the same universality class as the statistical model counterpart of the corresponding CFT. 
These height models have also loop representations. This similarity can be useful to get an 
idea about the loop properties of the statistical models with well-known minimal CFTs. 

4 Some Examples 

In this section we apply the method introduced in section 3 to the minimal conformal field 
theories with well defined fusion structure and also WZW SUk(2) models. We will also point 
on the consistency of these loop models with the Cardy's boundary states. These consistency 
is a hint to believe that it may be possible to extend the results in to the level of the boundary 
partition function [7J. For notational convenience in this section of the paper we will drop 
the hat of boundary operators. 

Ising model: The simplest example is the Ising model (^2,^43) since the model has 
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diagonal modular invariant partition function the fusion matrices of the boundary operators 
is the same as the bulk case. The fusion graphs are as fig 2 so the boundary states are as 
follows 



These equations reflect the Z2 symmetry corresponding to changing the sign of spin, this is 
also evident in the loop representation; n £ = n\ = 1. Both operators give k = 3, these loops 
are the domain walls between different spins. It is worth mentioning that this symmetry 
comes from the natural Z2 symmetry of Dynkin diagram. The operator a with n a = y/2 



describe the domain walls of Fortuin-Kasteleyn (FK) clusters. In the above calculation we 
considered only the largest eigenvalue of the fusion graphs, however, it is also possible to 
consider other eigenvalues as the weight of the loops, the cost is accepting complex local 
Boltzmann weights for the corresponding height model. Since loop models are generically 
non-local theories accepting complex Boltzmann weights is equal to accepting non-unitary 
theories. By this introduction one can accept the possibility of loop models with n = ±\/2, 
for the loop model corresponding to the A3 diagram of spin operator. 

Tri-critical Ising model: The next simple example is the tri-critical Ising model, 
(A3, A4) which we have diagonal modular invariant partition function. The boundary CFT 
of this model was discussed in |18j . There are 6 boundary operators 1, e, e' , e" , a and a' with 
the fusion graphs as fig 3 and the following Cardy states 



|1> = C[|l))+r ? |6))+r ? |6')) + |6")) + ^2k / )) + v / 2k))]; 

|e) = C[r ? 2 |l))-^ 1 | e ))-r ? - 1 |e'))+r ? 2 |e // ))-^2 ?? V / )) + v^ ?? -V))]; 

|e') = C[r ? 2 |l))-^ 1 | e ))-r ? - 1 |e'))+r ? 2 | e ")) + ^V / ))-^-V))]; 

e") = C [\l))+r ] \e))+r ] \e')) + \e''))-^W))-^2\cj))) ] 

a') = ^2C[|l))-r ? | e ))+r ? | e '))-k // ))]; 

\a) = ^C^ll^ + ^I^-^V))-^"))], (4-5) 




(4.4) 



corresponds to free boundary condition. The loops in the dense phase have k 
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where 




The boundary states corresponding to boundary 
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1 a' £" 
e o e' 

Figure 3: Graphs of fusion matrices of the boundary primary operators in the tri-critical 
Ising model in the upper row from the left to the right the fusion graphs of 1, e and e', in the 
lower row from the left to the right the fusion graphs of e" , a and a' . The fusion graph of e 
is A4 plus T2 , they are connected to each other by folding duality. The fusion graph of a is 
T 2 ® A 3 . 




operators 1 and e" can be transformed to each other by just changing the sign of spin 

operators, i.e. Z2 symmetry. They have also the same loop weight n = 1 comes from the 

largest eigenvalue of the fusion matrij^. The boundary states e and e' are connected also by 

just changing the the sign of spin states. The weight of the loops is n = 2cos(^) with k = 5 

in the dense phase. This loop model corresponds to the boundary of geometric clusters at the 

geometric critical point of tri-critical Ising model or Blume-Capel model |20| . The operator 

a' describes a loop model with n = y/2 corresponding to k = ^ in the dense phase which 

is related to the boundary of spin clusters and also vacancy clusters in Blume-Capel model 

|20t I21j. The interesting point for tri-critical models is the equality of critical exponents for 

spin clusters and FK clusters [20\ I21j. The operator a is related to the degenerate boundary 

condition and the corresponding loop model with n = 2\/2cos(^) is non-critical, however, 

it is easy to see that the fusion matrix of this operator is decomposable to simple matrices 

2 To get the loop weights we consider one simply connected part of the fusion graph as an adjacency graph, 
the other parts of the graph have always equal largest eigenvalues. One can see that these different parts are 
folding or orbifold dual of each other, see [19] 
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N a = Nt 2 <8> Na 3 so one can define for this graph two-flavor loop model with weights n\ = \[2 
and 77-2 = 2cos(^). One can conclude from the above discussion that those operators with 
the same loop representations are connected to each other by folding and orbifold duality 
and it is also possible to see these symmetries in the level of boundary states. 

Similar to the previous subsection one can also consider other possible loop weights come 
from the other eigenvalues of the fusion matrix. The eigenvalues of the fusion matrix of 
the operator e are n = ±2 cos (7^), ± 2 cos(-) anc ^ the eigenvalues of the fusion matrix of the 
operator a' are ±\/2, 0. The eigenvalues of the other operators are a subset of the above 
eigenvalues. Interestingly apart from the negative eigenvalues the above weights can be fitted 
with the boundary loop weights in [51 [6] . 

Three states Potts model: The next example is the first non-diagonal case, 3-state 
Potts model {A±,D±) with 8 boundary operators 1, ip, ip* , e, a, <r , and (j>2,2, see [T2"l 
[T3l I22j . The fusion graphs are given in fig 4. Following Cardy's argument one can show 
that the operators l,^,^ correspond to fix boundary conditions and the corresponding 
boundary states can be transformed to each other by Z3 symmetry, i.e. the symmetry of 
Dynkin diagram D4. They also have the same quantum dimensions n\ = = = 1. 
The operators e, a, describe the fluctuating boundary conditions [23J and all have the 
same kinds of fusion graphs with n e = n a = n a \ = 2cos(^). In the dilute phase one can 
consider the property of the curve. In the lattice 3-state Potts model these loops 

are the same as the domain walls of spin clusters. The fusion graph of the operator <j)\ 2 is 
two D4 graphs. This operator describes fix boundary condition and has loop model with 
n = y/3 which is equal to the loop model of domain walls in FK clusters of 3-state Potts 
model. The operator ^2,2 describes degenerate boundary condition and the corresponding 
loop model with n§ 2 2 = ^/ 9 +^v / 5 j g non-critical, however, decomposition is possible. In this 
case one can write Nfa 2 = Nt 2 <8> -^D 4 and so the corresponding two-flavor loop model has 
weights n\ = y/3 and rti = 2cos(|). 

The fusion matrix of e as was discussed in the case of tri-critical Ising model has the 
eigenvalues n = ±2cos(|), ± 2 cos(i) an< ^ * ne eigenvalues of the Njj 4 are ±v3, 0. These loop 
weights can be fitted with the boundary loop weights in [SJ [6]. 

Tri-critical three states Potts model: The next interesting example is tri-critical 3- 
state Potts model (-D4, A§) it has non-diagonal modular invariant partition function and also 
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Figure 4: Graphs of fusion matrices of primary operators in three states Potts model. In the 
upper row from the left to the right the fusion graphs of 1 and e, in the middle row from the 
left to the right the fusion graphs of 0i 2 and ^2,2 and in the lowest row the fusion graphs 
of tp and a. The fusion graphs of tp' and cr' can be derived from the fusion graphs of ip and 
a by the following exchanges ip ip' and a <-> o^. The fusion graph of e is plus two T2 
graphs, they are connected to each other by folding duality The fusion graph of ^2 is two 
D4 and the fusion graph of $2,2 is ?2 <8> D 4 . 
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it is not part of Pasquirer's A-D-E models. The boundary states of this model have not 
been investigated systematically so far. The boundary operators of this model are: 0j with 
i = (r,a), r = 1,2,3 and a = 1, ...,4. The fusion graphs for the boundary operators in this 
case are given in the Appendix. The boundary states corresponding to boundary operators 
01,1) 01,3; 01,4 can be transformed to each other by Z% symmetry of spin operators and should 
correspond to fix boundary conditions with n = 1. The operators 02,1,02,3,02,4 have also 
the same property with the same fusion graphs with n = 2cos(y). In the lattice tri-critical 
3-state Potts model they are domain walls of geometric clusters of geometric critical point 
|20j with k = 4g. The operators 03. i, 03,3, 03,4 can be transformed to each other again by Z% 
symmetry but they have loop weights bigger than two; n = 2.246. The operators 02,2 and 
03,2 have also loop weights bigger than two and related to degenerate boundary conditions. 
Finally the graph of 0i,2 is equal to three -D4 graphs with n = \^3. In the dilute phase this 
weight describes the domain walls of spin clusters in the lattice tri-critical 3-state Potts 
model with k = 4=. 

The fusion graph of 02, 1 is the sum of two graphs Aq and T3. The fusion matrix has the 
eigenvalues n = 2 cos(^) with j = 2, 3, 4, 5, 6. The eigenvalues of the fusion matrix of 0i,2 are 
n = ±\/3, 0. Interestingly again apart from the negative eigenvalues the above weights can 
be fitted with the boundary loop weights in The fusion graph of 02,2 is decomposable 

as T3 (g) D4 and so it is possible to define two crossing loop models in this case. The fusion 
graphs of 03, 1 is not decomposable to simple graphs so it is not possible to extract critical 
loops also for 03,3 and 03,4 which are in the same sector. Although the loops, extracted by 
our method, corresponding to the above operators are not critical but by considering the 
fusion graph of the ground state of the above adjacency graph it is possible to extract critical 
loops, we will not discuss this method here, for more detail one can see |16j . The fusion 
graph of 03,2 is decomposable but not to the simple graphs, i.e. A^ 3 2 = A T 2 (8>ATd 4 . Another 
possibility to extract critical loops for ^3,1 is by considering other eigenvalues of the fusion 
matrix of this operator. The eigenvalues of Nj r , , are ± sm ^l ) , ± sm \ I ) and ± s '°fxl , the last 

L b <P3,1 sm(-) ' sm(-) sm (iZL) ' 

two cases have critical loops. 

Minimal models: Finding loop models by the above method is completely general and 
applicable for more general cases. Take a pair {A, G) from the equation (|2.ip then it is pos- 
sible to correspond at least two different kinds of loop models for these minimal models with 
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the following weights 

n = 2cos(-), n = 2cos(^). (4.6) 

g h 

They are the largest eigenvalues of the fusion matrices of 0i 2 and 4>2,i- One can also consider 
the following SLE drifts for these loop models 

K = 4^, K = 4^A_, g -h = -1. (4.7) 

The other eigenvalues of G can be written as 

n = 2cos(— ), (4.8) 
9 

where m is one of the Coexter exponents of the graph G. they are listed in the table 1. 

It is possible to consider loop models for the above eigenvalues as before, however, they 
are not still all the possible loop models because as we already showed in some cases one can 
define two flavor loop models for decomposable fusion graphs. It is also possible as the case 
of the fusion graph of 03 5 i in tri-critical 3-state Potts model to have matrices with relevant 
non-largest eigenvalues. We believe that they are relevant because the same loop weights 
appear in the classification of Jacobsen and Saleur [5]. 

Although so far we have given more familiar examples as the possible candidates for our 
loop models but it is also possible to extract systematic examples for the above proposals 
by using Pasquier's ADE models and Dilute ADE models [10[ 111]. Pasquier's ADE models 
give a lattice realization for the (A, G) series with g — h = 1 and the description briefly is as 
follows: define an RSOS model by using the graph G this height model at the critical point 
can be described by a the minimal CFT then map this height model to loop model [21] at 
the critical point with n = 2cos(^) which is the same as the loop model that we proposed 
in (|4.6p . Of course the method proposed in this article and [3] is highly influenced with 
Pasquier's ADE models but it has something more to say by connecting the loop properties 
to the fusion properties of the primary operators. To get the dilute loop models and the 
loop models corresponding to tri-critical models we need to use Dilute ADE models. These 
models have rich phase diagrams with four branches: branch 1 and 2 have central charges 
c = 1 — g(gj-i) an d branch 3 and 4 have c = | — ■ One can also map this height models 
to 0(n) loop models with the non-intersecting bonds on the square lattice with the partition 
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1 u 




Figure 5: The Boltzmann weights of the different vertices in the 0(n) model on the square 
lattice. 



function 

Z = Y / u Nu v N ^w N -n N , (4.9) 

where the weights for different plaquettes are given in Fig 5 and the N U ,N V and N w are 
the numbers of different plaquettes [25]. This generalized 0(n) loop model apart from the 
critical properties at u = w = \ an d v = has four other branches coincide with the four 
branches of dilute ADE models [27]. The weights are given by 

n = -2cos(26>), 

1 

111 = 

2 - [1 - 2 sin(0/2)] [1 + 2 sin(#/2)] 2 ' 
u = ±4u7sin(6»/2)cos(vr/4 - 0/4), (4.10) 

v = ±u;[l + 2sin(0/2)], 

where ^ < 8 < it, < 9 < ^, — ^ < < and — ir < 9 < — ^ are the intervals corresponding 
to branches 1, 2, 3 and 4 respectively. They coincide with the different branches in the 
dilute ADE models. 

It is interesting to investigate the connection of the above loop model to the SLE. There 
are different methods to do that here we use the magnetic operator to find the SLE drift. It 
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was shown in [25] by the numerical calculation that the magnetic exponent of the branch 1 



and 2 is identified with 2h where 

2 ' 2 



_ ((m + l)r - ms) 2 - I 
4m(m + 1) 

and m is related to the central charge of the theory by c = 1 — m ^ • Its connection to the 
loop variables comes from the relation ^ + ^ — 2 = m( . derived from the coulomb gas 
method [25]. The connection of the magnetic exponent to the SLE drift is as follows 



2h„. n „„ = < 8 -"» 3K - 8 >. (4.12) 
-3->— 32k v 1 



Using the above equation the SLE drift at the branches 1 and 2 of the loop model (|4.9p can 
be derived as follows 

K = -• (4.13) 

7T 

This result is consistent also with our expectation from the second level null vector of minimal 
models [28], it is also consistent with the recent direct investigation by using holomorphic 
variables |29] . 

Back to the height model representation one can summarize following results: the branch 
2 of the ADE models corresponds to the dilute loops with n = 2cos(^) and the branch 1 is 
the dense phase of tri-critical models with n = 2cos(^). The results for some of the simple 
cases are as follows: 



branch 2: 


A 2 


= critical percolation, 


c 


= 


n 


= i; 


branch 1: 


A 2 


= critical Ising 


c 


= 1/2 


n 


= i; 


branch 2: 


A 3 


= critical Ising 


c 


= 1/2 


n 


= V2; 


branch 1: 


A 3 


= tri-critical Ising 


c 


= 7/10 


n 


= V2; 


branch 2: 


A, 


= tri-critical Ising 


c 


= 7/10 


n 


= 2 cos (-) 
5 


branch 2: 


D 4 


= critical 3-state Potts 


C : 


= 4/5 


n - 


= v^; 


branch 1: 


D 4 


= tri-critical 3-state Potts 


c 


= 6/7 


n 





Using the above method it is easy to find the lattice realization for most of the proposed 
loop models, the results are interestingly consistent. Following the same method it is possible 
to extract the loop models corresponding to minimal CFTs, however, the loop model for the 
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non-diagonal cases with g — h = — 1 is not extractable with this method because we are not 
able to find the dense phase of loop models for these cases. It seems that the dense lattice 
height model has not been proposed for this case. 

To conclude this subsection we proposed some loop representations for the minimal CFTs 
by using fusion of boundary operators. Then since ADE models give a lattice statistical 
model representation for minimal CFTs we used these models to extract physical loop mod- 
els corresponding to ADE models. The fractal properties of these lattice loop models are the 
same as the loop models that we proposed by using the fusion of primary operators. 



SUk(2) Models: It is possible to follow the same calculation for every unitary minimal 
model. For example for WZW SUk(2) models the classification of modular invariant parti- 
tion functions is based on A-D-E-T graphs with g = k + 2. The same method as the minimal 
models is applicable here and one can find boundary operators <pj with 1 < j < k + 1. The 

sin( ) 

loop models have weights dj = sin (I) ■ Only j = \ has critical loop representation with the 
following loop weight 

n = 2cos(-^), (4.14) 

with k = 4|±| and k = 4|±| for the dilute and dense phase respectively. The other loop 
models are not critical except for k = 4 with n = 2. The fusion graphs of the operators 
with j | is not decomposable to the simple graphs, however, the non-largest eigenvalues 
can be still relevant. For example take k = 5 with j = 3/2, the fusion graph is similar to 
the one part of the (f>si fusion graph of the tri-critical 3-state Potts model, the right one in 
fig 6. The eigenvalues are ± sm \ I } , ± sm ^ I ) and ± Sm /xl , the last two cases have critical 

to to sm(-) ' sin(-) sm (4ZL)' 

loop representation. The similarities between fusion graphs of SUk(2) models with minimal 
models is not just an accident they are based on the coset construction of the minimal models. 



5 Discussion 



We proposed a method to classify some possible loop models consistent with the conformal 
boundary conditions for generic rational CFT: take the simply laced classification of the 
corresponding minimal CFT then find the boundary operators and also the fusion matrices, 
make the 0(n) loop model of the primary operator by the method that we discussed in section 



17 



3 and [4 J . We think that there should be some connections between these loop models and the 
SLE interpretation of CFT investigated in [28J which is based on the connection of SLE with 
the null vectors in the CFT. This connection is not complete even for minimal CFTs because 
we do not know how to explain the boundary operators with the same loop model but with 
the different null vectors, for example in the three states Potts model e, a and at are in the 
same sector from the boundary CFT point of view but just e and a have the required second 
level null vectors. However, from null vector point of view this correspondence is not clear 
but it is possible to show that in the partition function level this similarity is more known. 
Another way to look at the results of this paper is by conjecturing the largest eigenvalue of 
the fusion graph as the possible loop weight for the loop model in the universality class of 
the corresponding CFT without defining any height model on the fusion graph. 

One possible generalization of the above construction is by considering graphs with largest 
eigenvalue bigger than 2 as an adjacency graph of fused RSOS model and then extracting the 
loop model by the method investigated in [16]. The other interesting direction is to investi- 
gate the modular invariant partition functions of loop models and their possible connections 
to the classified modular invariant partition functions of minimal models, this is related to 
investigate more directly the connection of our method to the classification of [5, 6j. 
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6 Appendix 

In this appendix we list the fusion graphs of the boundary operators in tri-critical 3-state 
Potts model. The fusion graphs are given in fig 6. The fusion graph of 0i 5 4 can be derived 
from the fusion graph of the operator 0i 5 3 by the following transformations: 

01,3 <-> 01,4 02,3 <-> 02,4 03,4 <-> 03,3- (6.15) 

The fusion graph of 02,3 can be derived from the fusion graph of the operator cp2,i by the 
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...... 

#1,1 #1,2 fb #,4 #2,1 #2,: 

...... 

#2,3 #2,4 #3,1 #3,2 #33 #3/ 



#1,1 #1,2 #,3 #2,1 #2 

#3,3* 

#3,1 #3,2 #3,4 



G 



. . . . . . 

#1,2 #,4 #2,2 #2,3 #3 2 #3,3 



#1,1 #,3 & 4 



?3,1 P3.4 



4. 



■8- 



#,3 #2,4 #3,3 #3, 



G 



fa? 
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following transformations: 



01,3 01,1 02,3^02,1 03,3^03,1- (6.16) 

Finally the fusion graph of 02,4 can be derived from the fusion graph of the operator 02 i by 
the following transformations: 



'1,4 9>1,1 92,4 <-» <f2,l <P3,4 <P3,1- 



(6.17) 



To get the fusion graphs of 03,3 and 03,4 from the fusion graph of 03, i one just need to use 
the transformations (|6.16p and (|6.17p respectively. 

We shall call the part of the fusion graph of 03, i with two neighbor blobs T| , the lower 
index is the number of nodes and the upper index is the number of blobs attached to the 
neighboring nodes of the graphs starting from one of the extremes. These kinds of fusion 
graphs appear also in the fusion graph of 4>j=i of SU2(k) models. 
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